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A concise, somewhat personal, review of the problem of superfluidity and quantum 
criticality in regular and disordered interacting Bose systems is given, concentrating 
on general features and important symmetries that are exhibited in different parts of 
the phase diagram, and that govern the different possible types of critical behavior. A 
number of exact results for various insulating phase boundaries, which may be used to 
constrain the results of numerical simulations, can be derived using large rare region 
type arguments. The nature of the insulator-superfluid transition is explored through 
general scaling arguments, exact model calculations in one dimension, numerical results 
in two dimensions, and approximate renormalization group results in higher dimensions. 
Experiments on *He adsorbed in porous Vycor glass, on thin film superconductors, and 
magnetically trapped atomic vapors in a periodic optical potential, are used to illustrate 
many of the concepts. 

Keywords: Boson superfluidity; disordered systems; quantum phase transitions. 



1. Introduction 

It has been approximately 20 years now since the flowering of interest in the dirty 
boson problem. The problem is deflned, generally, as the nature of the insulating 
and conducting phases, and the phase transitions between them, in a system of 
interacting bosons in a random potential at zero temperature. A phase transition 
occurring at zero temperature, as a function of some auxiliary control parameter, 
such as density or magnetic field, is known as a quantum phase transition (QPT) 
since the particle dynamics are provided purely by the fluctuations in the ground 
state wavefunctionjflln a Bose system the conducting phase is believed always to be 
a superfluid (SF) (see, e.g., Ref. [4|), and this strongly distinguishes the dirty boson 



''At finite temperature, through standard arguments (see, e.g., Ref. |Tj) quantum mechanics is 
irrelevant near criticality, and the physics reduces to that of the corresponding classical model — 
the XY-model in the case of superfluid *He. The effects of disorder are then governed by the well 
known Harris criterion (Ref. |2}. For discussion specific to bosons on a disordered substrate, and 
further references, see Ref. \3\ 
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problem from the correspond rng A nderson localization and metal-insulator transi- 
tion problems in Fermi systemsl^^In particular, the conducting phase has a natural 
order parameter, absent in an ordinary metal, associated with the off-diagonal long 
range order of the superfluid. Moreover, in the absence of Pauli exclusion, repul- 
sive interactions are essential in forestalling condensation of a macroscopic number 
of particles into the single lowest localized free particle eigenstate of the random 
potential. There is therefore no sensible noninteracting limit about which to per- 
turb to study the physical finite density phases. Thus, although there are some 
analogies between the Fermi and Bose phenomena, the underlying physics is very 
different, and the two problems require entirely different theoretical approaches to 
their solution. 



1.1. Some "ancient" history 

Perhaps the earliest relevant experiments, which certainly were key to motivating 
early work on the dirty boson problem,^^ were on superfluidity of very thin films of 
■*He adsorbed in porous Vycor glass.^^ Vycor glass has a highly connected, 40% open 
(well above the percolation threshold), 3D sponge-like structure with pore sizes in 
the 4-8 nm range. In an interesting historical reversal, the experiments themselves 
were motivated not by disorder effects, but by the claimed earliest experimental 
evidence for an essentially ideal Bose gas, a do zen years in advance of its observation 
in magnetically trapped atomic vaporsP^^^ Of course, the Vycor system is not 
nearly as "clean" as the atomic system (in senses that should become clear below) , 
so the claim itself is not as clean as might be desired]^ 

Some of the Vycor superfluid density (ps) vs. temperature datePare reproduced 
in the upper left panel of Fig. [1] The ideal Bose glass claim emerged from an 
examination of the evolution of the shapes of the Ps{T] p) profiles in the vicinity of 
the transition Tc{p) as the overall helium density p was reduced. It was observeci^ 
that the exponent v describing the superfluid onset power law, ps ^ {tl"" below Tc, 
where t ^ {T — Tc)/Tc is the reduced temperature, showed a clear crossover from 
the classic bulk value, v ~ 2/3, observed at higher coverages, toward the ideal Bose 
gas value, vq — 1, at the lowest coverages. Quantitative theoretical predictions 
for th is crossov er, based entirely on a clean, non-disordered "effective medium" 
model, P l ^^ l ^'^ l were generally consistent with this scenario. 

Although the trends in the vicinity of Tc support an ideal gas interpretation, 
other elements of the data definitely do not, and this will serve as motivation for 
the remainder of this article. Motivated by QPT ideas, the remaining panels of Fig. 
[T] show plots involving various zero temperature quantities.^l^J In the upper right 
panel is plotted the extrapolated zero temperature superfluid density ps{0] p) vs. p 
(the horizontal axis should be divided by the 0.89 cm'^ total glass-plus-pore volume 



''The essential role of Vycor is that it acts to screen the long-range attractive part of the helium- 
helium pair potential, which in the bulk causes the vapor to condense into a dense fluid, preempting 
observation of the dilute regime in pure, bulk ^He. 
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Fig. 1. Superfluidity of helium in Vycor. Upper left: Superfluid density (proportional to the 
torsional oscillator period shift AP due to the decoupling of the superfluid from the oscillating 
substrate) data for ^He adsorbed in porous Vycor glass. The different curves correspond to different 
fixed values of the overall filling. The trend with decreasing coverage p, towards a more linear, 
less steep, onset just below the transition temperature Tc{p), is evident to the eye. Upper right: 
Extrapolated zero temperature superfluid density vs. coverage, with fits to ideal Bose gas and 
QPT-motivated functional forms. Lower left: Extrapolated zero temperature superfluid density 
vs. Tc compared to ideal gas models. Lower right: QPT-motivated power-law fit to the zero 
temperature superfluid density vs. Tc data. 



if one is fussy about units). Zero temperature superfluidity evidently disappears 
for p < pc — 4.35 mmol, corresponding to a coverage of about 1.5 monolayers 
(over the estimated 156 m^ pore surface area). If one interprets pc to be an inert 
background upon which the superfluid component "floats," then an ideal gas model 
would predict a linear relation Ps{0', p) (x p — pc- Evidently, one may force a linear 
fit to the lowest density data by adjusting the slope using an effective mass free 
parameter. The fitted value m*/m ~ 3.5 is at least consistent with one's intuition 
that interaction with the substrate would tend to increase the mass. However, 
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a more convincing fit, over a much wider range, is obtained with a superlinear 
fit ps(0;p) cx (p — Pc)", which serves to define our first QPT critical exponent 
w~ 1.7 ±0.3. 

The bottom two panels of Fig. [1] compare the rates at which Ps{0; p) and Tc{p) 
vanish with p. An ideal gas model predicts Tc cx p''/^, with spatial dimensionality 
d — i here. The lower left panel shows that one may force such a proportionality 
over the lower coverage data, again using an effective mass free parameter. The fact 
that the inferred value m*/m ~ 1.5 is different from the previous one should lead 
to some discomfort. Once again, as seen in the log-log plot of the same data in the 
right panel, a much more satisfactory power-law fit Ps(0;p) oc TcipY is obtained, 



1.2. Outline 

With the helium in Vycor system serving as an introduction, the outline of the re- 
mainder of this paper is as follows. This review is not intended to be comprehensive, 
and the author apologizes in advance for his personal choices as to which develop- 
ments deserve emphasis. In Sec. [2] the theoretical approach to the problem is intro- 
duced through the Bose-Hubbard and Josephson junction lattice models. In Sec. [3l 
the phases, phase diagrams, and accompanying phase transitions are motivated and 
summarized. In Sec.|4l exact results for the incompressible phase boundaries, based 
on the existence of exponentially rare, but arbitrarily large, disorder free regions, 
are derived. In Sec. [5] the quantum model is mapped onto a classical field theory 
Lagrangian. The various parameters appearing in the latter allow one to transpar- 
ently identify the different possible symmetries of the model, and accompanying 
QPT universality classes. In Sec. [H QPT critical scaling phenomenology is intro- 
duced, which, in particular, allows one to express the exponents u and 9 in terms 
of more familiar ones. Scaling also predicts the existence of universal sheet conduc- 
tances in d = 2, which has some experimental support. Renormalization group and 
epsilon expansion approaches to detailed model calculations are discussed in Sec. 
[71 The latter are poorly controlled compared to their classical counterparts, but 
nevertheless provide an compelling global picture of the stability, bifurcation, and 
merging of critical fixed points as a function of dimension that is consistent with 
the phenomenological scaling arguments. The review is summarized and concluded 

'^Why, over a large range of coverages, is the behavior near Tc consistent with an ideal Bose gas 
crossover, while over the same range the zero temperature quantities show strong deviations? The 
answer is not entirely clear but, firstly, critical scaling near Tc is entirely separate from QPT scaling 
near T = — for a detailed discussion of this in the context of the clean system, see Ref. il4l Thus, 
disorder is expected to have distinctly different impacts on the two ends of the superfluid density 
profile. Second, as discussed in Ref. jS] its manufacturing process makes Vycor very uniform on 
multi-pore scales, which aids the uniform effective medium approximation which forms the basis 
of the ideal gas model. However, one eventually expects to see strong deviations, as the disorder 
begins to impact the ideal Bose gas physics, presumably at yet smaller values of p — pc that are 
beyond experimental resolution. 
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in Sec. E 

2. Bose-Hubbard and Josephson junction Models 

In analogy to the Hubbard model of electronic propagation in a crystalline solid, a 
useful starting point for the study of bosons in a disordered medium is the Bose- 
Hubbard model of lattice bosons. The Haniiltonian is 

Hb ^ ~J^^ - ~ /"X/"* + ^^^^f>'ii''^i - ^) + Uifii, (1) 

i i i 

in which aj , di are boson creation and annihilation operators on site i (taken to lie on 
a regular lattice, and obeying the usual Bose commutation relations [a^, Oj] ~ Sij), 
fii — aldi is the (non- negative integer) site number operator, the first sum is over 
nearest neighbors, J — /2m*<P is the hopping amplitude, where d is a measure 
of the pore diameter and m* is a boson effective mass, is the chemical potential, 
> the on-site repulsion, and Ui the random site energies. In this model each site 
represents a pore, and the microscopic structure of the individual pores has been 
subsumed into the model effective parameters. The same picture obtains for bosons 
hopping between minima of the opti cal latt ice potentials that are used to perturb 
magnetically trapped atomic vapors.^I^H^ Clearly, a more realistic model would 
include disorder in J, V ^ and the lattice site positions, but, as in the fermion case, 
these add nothing new to the basic physics implied by the random ui. One could also 
add higher order Hubbard bands, corresponding to higher order intra-pore single 
particle excited states, but at low temperatures, lacking Pauli exclusion, these will 
not be populated, and can be ignored. The precise probability distribution of the 
Ui is not important, but it is useful to consider them as independent from site to 
site, with an even distribution (in particular, with mean zero, which normalizes the 
chemical potential), and with support on a finite interval — A < Ui < A, with the 
bound A serving as a free parameter. 

A closely related model is the Josephson junction, or quantum rotor model, 

TYj = - ^ Jy cos(0i - (f)j) - fJ-'^fii + ■^V'^hf + '^u^hi, (2) 

(ij) i i i 

in which hi, (pi are conjugate number and phase operators, defined by the commu- 
tation relations [(pi, hj] — iSij. In this model hi can take negative as well as positive 
integer values. The mapping between the two models is provided by identifying 
aj h^^e^'^^di ^ e^'"^^hy^. This Hamiltonian originated as a model of granular 
superconductors, with hi representing the number of Cooper pairs on grain i, and 
with Josephson couplings Jy between the phases of the neighboring grains. The 
hopping term clearly lowers the energy when neighboring phases are aligned, and 
the superfluid, or superconducting, state corresponds to the appearance of long 
range order, and corresponding nonzero local averages (e**^' ) ^ [corresponding to 
the usual Bose condensate (di) ^ in Quantitative agreement between the 
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two models is obtained in the limit of large mean site occupancy rii — (hi) 3> 1, 
where number fluctuations can be neglected in the hopping term and one identifies 

By separating amplitude and phase in this way the underlying symmetries of 
the model, which are critical to identifying the possible QPT universality classes, 
are much more transparently exhibited. First note that under an integer transla- 
tion fii fii + no, which preserves the commutation relations, the Hamiltonian 
transforms in the form, 

H,7(m) ^ H.7(/i - noV) + (^-/ino + ^^^"o^ Nl, (3) 

where Nl is the number of lattice sites. It follows that the phase diagram is periodic 
in n with period V (see Fig. 01 the additive constant serves simply to increase the 
overall density by hq). Second, the particle- hole transformation corresponds to the 
reflections, 

-n,, (f)i (4) 

which also preserves the commutation relations. Under this transformation, 
7i(yLt, {ui}) T-L(—^, {— lii}). Together with the translation symmetry, this implies 
that at the special values ^ = noT^/2, in the absence of site disorder, Wi = (but ar- 
bitrary Jy ), the Hamiltonian is invariant (up to an additive constant), implying the 
existence of a special symmetry between particle and hole excitations. Furthermore, 
for symmetrically distributed Ui, the transformed Hamiltonian is statistically iden- 
tical to the original. In either case, the phase diagram has an additional reflection 
symmetry about the points ^ = no 1^/2. 

These special particle-hole symmetric and statistically particle-hole symmetric 
points (which T-Lb can only exhibit indirectly — see below) will play an important 
role in understanding the critical behavior. Furthermore, the fact that the hopping 
is unaffected by the transformation motivates consideration of disorder in both 
Jij = J{\ + SJij) (which serves to define an overall control parameter J) and 
Ui in this model. Understanding the effects of breaking these symmetries (which 
ultimately leads back to Hb, but now understood in a broader context) will turn 
out to be key to understanding the global renornialization group flows (SeclTj). 

3. Phases and phase diagrams 

By considering the energy balance between the various terms in ([T]) and ([2]) one can 
understand the basic topology of the phase diagrams. First the results for the Bose- 
Hubbard model (H]) will be motivated and presented. Additional features special to 
the Josephson junction model ^ will then be discussed. Detailed calculations can 
be found in the literature: the original mean field predictions for the phase diagram 
were made in Ref. [7J subsequent perturbation calculations can be found in Ref. [531 
[231 and numerical simulations can be found, for example, in Refs. [T71 [IHl [IHl [23 

[mm 
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Fig. 2. Left: Schematic illustration of lattice bosons near unit filling n = 1 at some value of the 
hopping amplitude in the interval < J < Jc- In the Mott insulating phase (n = 1), the effective 
wavefunction of each particle spreads only a finite distance ^(J), and the state is insulating. For 
n > 1 (n < 1), the extra particles (holes) travel freely within the essentially inert background, and 
the state is superfiuid for arbitrarily small \n — 1|. Right: Associated phase diagram with Mott 
insulating (MI) and superfiuid (SF) phases. The MI are incompressible, with fixed integer filling 
over a finite range of /i. The SF is compressible, with integer filling only along some set of lines 
extending from the Mott lobe tips. The phase transition through the tip, where a "hidden" particle- 
hole symmetry is asymptotically restored, is in the universality class of the {d + l)-dimensional 
classical XY model. Elsewhere, the transition is equivalent to onset of superfiuidity at zero density 
in a continuum dilute Bose gas. In this sense, as indicated, "periodic Vycor" would indeed produce 
the desired ideal Bose gas crossover upon approach to integer filling. 

3.1. Clean system 

We begin with the clean, non-disordered system, Ui = 0, corresponding to bosons 
moving in a perfectly periodic background potential. The left panel of Fig.[2]sketches 
the essential physics. If there is no hopping, J = 0, then the sites decouple, the 
site occupancies are good quantum numbers, and there are exactly uq particles on 
each site for no < ^i/V < no + 1. There is a unit jump in filling at each integer 
value of pl/V . For finite J, particles are no longer confined to a single site, but 
may wander to neighboring sites. However, for small J /V the mutual repulsion 
produces an energy gap of order V whenever two particles occupy the same site. 
So, although virtual exchanges do occur, the net effective distance ^(J, no) a given 
particle wanders remains finite. There is a corresponding energy gap /x±(J, no) = 
±[i?(J, noiVi ± 1) — i?(J, noiVi)], where E{J,N) is the ground state energy for N 
particles, for adding (-I-) or removing (— ) a particle. Thus, for chemical potential in 
the interval H-{J, no) < /i < fJ-+{J, no), the ground state wavefunction is insensitive 
to /i, and, in particular, the density remains fixed at no- This incompressible phase 
is known as a Mott insulator (MI), analogous to the corresponding electron half- 
filled band insulators, and the resulting lobe structure of the phase diagram is 
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sketched in right panel of Fig. [2l All of these notions can be confirmed through 
direct perturbation theory in j/^}2312l] 

Only for sufficiently large J is the hopping sufficiently vigorous to overcome the 
repulsion: at a critical value J Jc{no) the hopping range diverges, ^ oo, and 
the gap closes, e± 0. For J > Jc one enters the superfluid phase. One may define 
a quantum critical exponent correlation length exponent v via ^ ~ \J — Jc\~'^ , and 
dynamical exponent z via e± ^ | J — Jc\^'^ . It will be shown in Sec. [5] that this 
transition, at which, as alluded to above, an effective particle-hole symmetry is 
restored, is in the universality class of the classical {d + l)-dimensional XY model: 
v = h'^^i, z = 1. Thus, for example, in d = 2, vf^ ~ 0.671 should be the same as 
that measured at the ^He lambda transition. 

Between Mott lobes, /^+(J, ng) < /i < /i_(J, no + 1), the density varies contin- 
uously between the two integer values, uq < n{fi, J) < uq + 1. Near the Mott lobe 
boundaries one may think of the difference n — uq as extra particles (or riQ + l — n as 
extra holes) propagating atop the Mott phase background. Now lacking a barrier to 
hopping between sites, these particles (or holes) form a superfluid (lower left part of 
Fig. [5]) . Although the notion of coexisting insulating and superfluid components is 
misleading (bosons are identical particles, and there is strong exchange between the 
two groups), the effective theory near the phase boundary (analogous, perhaps, to 
the Fermi liquid theory of effectively free electron excitations near the Fermi surface) 
is indeed that of a dilute Bose gas. The exchange with the background insulator 
gives rise to a strongly renormalized (especially near Jc) effective single particle 
mass m* . In this sense, if Vycor had perfectly periodic pores, experiments would 
indeed see a crossover to dilute Bose gas behavior, with pc = no/d'^ corresponding 
to the density at integer filling. Perhaps the closest experimental re alizati on is that 
of trapped atomic vapors in an additional imposed optical lattice.l^ ^ l ^^ l However, 
although both MI and SF phases have been detected in these experiments, trap 
inhomogeneities, finite size effects, and other measurement constraints, are such 
that detailed critical phenomena are not yet resolvable. 

As stated, the very different critical behaviors at the "commensurate" transition 
at Jc, and the line of "incommensurate" transitions below Jc are indicative of a 
"hidden" particlc-hole symmetry at the former. The lack of an exact symmetry in 
Ti.B means that this symmetry is effectively restored at the nontrivial point where 
fi+{Jc, no) = fi^(Jc, hq) (analogous to the asymptotic restoration of the underlying 
Ising up-down symmetry of density fluctuations near a liquid- vapor critical point). 
During the transition to superfluidity at integer filling, particle and hole excitations 
exist in equal numbers, tunneling through each other, and achieve superfluidity 
simultaneously. Below Jc the transition takes place by feeding extra particles or 
holes into the system, explicitly breaking this symmetry, and only the predominant 
excitation achieves superfluidity. 
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n=0 

Fig. 3. Left: Schematic illustration of lattice bosons near unit filling in the presence of bounded 
site disorder. If the disorder is not too strong (5 = A/V < there is still a (shrunken) Mott 
lobe with a finite energy gap for adding or removing particles. Unlike in the pure case (Fig. 
[JJ, superfiuidity is not generated immediately outside this lobe. For sufficiently small n — 1, the 
additional particles are Anderson localized by the residual random background potential of the 
effectively inert layer. A finite compressibility distinguishes this new insulating Bose glass (BG) 
phase from the MI. The superfluid critical point fia{{J) occurs only once the added particles have 
sufficiently smoothed the background potential that its effective lowest lying single particle states 
become extended. Right: Associated phase diagram, with MI, BG and SF phases. The transition 
to superfluidity is always from the BG phase, is in the same universality class along the entire 
transition line, and is ultimately the correct description of helium in Vycor. 

3.2. Disordered system 

Consider now the addition of site disorder. The left panel of Fig. [3] motivates the 
existence of a third phase, the Bose glass (BG) phaseO that intervenes between the 
Mott and superfluid phases.'^ Beginning again with the J — limit, for sufficiently 
bounded disorder, A < V/2, there remains an interval (no — 1)V^+A < /i < noV~A 
over which every site still has exactly no particles. However, for /z just outside this 
interval, some sites with have an extra particle (or hole), and the fraction of such 
sites will vary continuously with /x — even at J = the state is now compressible. 

For small J, there will still be an interval fi~ (J, A, uq) < fi < ^i^{J, A, tiq) where 
mutual repulsion continues to the dominate the disorder, and the incompressible 
Mott lobe, though shrunken, still survives (right panel of Fig. [3]). A rare region 
argument will be used in Sec. |4]to show that /i±(J, A, ng) = /^±(^, 0, ng) T ^ are 

"^The term "Bose glass" may have originated in the early work of Ref. 1251 who studied a Hartree- 
Fock approximation to JTJ. It transpires that such an effective single-particle treatment of the 
interactions is a very poor approximation for disordered bosons, and their main conclusions were 
disputed in Ref. 1261 
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simple translates of those of the clean system (Fig. [5]). 

The compressible phase just outside the Mott lobe, however, can no longer 
be a superfluid. The extra particles (or holes) propagating atop the background 
Mott phase now encounter a background random potential, and the usual Ander- 
son argumenti^ show that the low energy effective single particle states must be 
localized, and the state therefore remains insulating. Only after a sufficient density 
of particles has been added, fi > ^sfiJ), is the residual random potential sufficiently 
smooth that that the lowest lying state becomes extended, and the transition to 
superfluidity takes place. Note the key role here of the pair interactions that allow 
particles to gradually "fill in" the deeper minima in the background potential. 

Note that for A > V/2 (which includes unbounded, e.g., Gaussian, disorder) the 
Mott lobes are entirely destroyed, and only the Bose glass and superfluid phases sur- 
vive. This is presumably the case for helium in Vycor, where the pore irregularities 
are too strong to permit a Mott phase. 




Fig. 4. Schematic phase diagrams for the Josephson junction model (adapted from Ref. \27l . 
Periodicity in fi/V is evident. Even site energy distributions have also been assumed, hence the 
symmetry under fi — > — /i as well. Left: Clean case. Center: Generic site-disordered case, with 
5 = A/V. Right: Case of pure hopping disorder. The compressible Bose glass (BG) is replaced an 
incompressible random rod glass (RRG) on the special lines at integer filling. The penetration of 
the superfluid phase all the way to J = at half-integer filling may be understood via a mapping to 
a sp in- ^ quantum planar XY model with random exchange, but vanishing out-of-plane magnetic 
field.EZl 
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3.3. Phase diagrams for the Josephson junction model 

For comparison, phase diagrams for the Josephson junction model are shown in 
Fig.m The clean and site disordered models (left and center panels) are essentially 
identical to their Bose-Hubbard counterparts (Figs. [2] and [H respectively), except 
for the periodicity in ^/V, and symmetry under inversion /i for evenly 

distributed Ui (reflecting the existence of an exact, rather than effective, particle- 
hole symmetry). 

On the other hand the case of pure hopping disorder (right panel of Fig. 2]) 
has no Bose-Hubbard analogue. In the absence of an exact particle-hole symmetry, 
random Jy always produces, in a renormalization group sense, an effective ran- 
dom site energy even if none is present in the original model — hence the choice of 
nonrandom J in ([T]). The exact particle-hole symmetry for integer and half-integer 
/x/y drastically changes the character of the model along these lines. For integer 
fi/V — no, although the Mott gap closes at a finite value J = Jm, the random 
rod glass (RRG) phase (see Sec. [5]) beyond it is incompressible. However, instead 
of vanishing over a finite interval, for J > Jm the compressibility is finite for any 
fj,/V 7^ no (thus indicating a Bose glass phase}, and vanishes with an essential sin- 
gularity, ~ e-'=o/\tJ-/v~no\ ^ \^J■/V — no\ — > 0.123 There are corresponding stretched 
exponential tails, as oppo sed to the Mott phase linear exponential tail, in the tem- 
poral correlation function.123 

At half- integer ii/V the particle- hole symmetry has a different effect, eliminating 
the glassy phase for all nonzero J. The model may be mapped to a spin-i quantum 
XY model, with the the Ui corresponding to mean-zero z-axis magnetic fields. If 
Ui ^ 0, then for sufficiently small J it is energetically favorable for a large fraction 
of the spins to align with the field along z, destroying the in-plane ferromagnetic 
order that signifies superfiuidity — this is the Bose glass phase. Only if = can the 
in-plane order, at least in higher dimensions d >1@ be maintained for arbitrarily 
small exchange coupling. 

4. Droplet model of the Bose glass phase and the MI BG phase 
boundary 

Before describing more formal theoretical approaches, I turn here to some very 
simple arguments that provide a more intuitive picture of the Bose glass phase, 
determine the exact vertical width of MI-BG phase boundary (Fig. [5]), and establish 
essentially rigorously that there can never be a direct MI-SF transition. As in the 
clean case, for given J, hq the Mott phase is insensitive to fj., and described a unique 
ground state wavefunction. The phase boundary again occurs when /i coincides with 
with the lowest lying single particle (or hole) excitation ^±( J, A, no) = ±[iJ(noiVL± 
1, A, J) — E{nQNL, A, J)]. The nature of these excitations, and thereby the values 

"For a strong-disorder approach to the ID problem, complementary to the weak disorder limit 
treated by Kosterlitz-Thouless-type theories, see Refs. [28] 1291 
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Fig. 5. Sketch of the Mott lobe boundaries as a function of the disorder (adapted from Ref. \27l . 
For simplicity, we consider here the ng = lobe of the Josephson junction model, for which 
= — /i— . The phase boundaries are simply those of the pure system translated inwards by 
the disorder distribution half- width A: fi±{J, A) = fi±{J,0) =F A. The result, for example, is a 
corner, or slope discontinuity, at the Mott lobe tip Jjv/(A) satisfying /i-|_(Jjv/, 0) — fi—(JM , 0) = 2 A 
(marked by stars). 



of fi±, will now be established. 



4.1. Low-lying particle and hole excitations 

Consider the addition of a single particle to the system. The energy required can 
never be less than it would be if all site energies were as small as possible, ui = — A 
for all i. This situation is equivalent to a simple shift, fi fj, — A, and the bound 
/Lt_|_( J, A, no) > /i+(J, 0,no) — A follows immediately. On the other hand, with an 
exponentially small density, scaling as p{R) ~ e"'^/^"^''-" , where Ro is some length 
scale depending on the precise distribution of the Ui, there will be rare, isolated 
"droplets" enclosing a sphere of arbitrarily large radius R, on which all Ui < — A-|-i5w 
are nearly uniform, and within an arbitrarily small energy Su of —A. The single 
particle excitation energy of such a droplet (above the uniform Ui — —A value) 
scales as Se ~ ciSu + C2J {d/ R)^ , where ci, C2 are appropriate geometrical constants, 
and the second term represents the ground state of a particle in a box of radius 
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In an infinite system, both terms can be made as small as one wishes, and 
one obtains the opposite inequality /i_|-(J, A, ng) < /i_(-(J, 0, ng) — A. One therefore 
concludes that ^+(J, A,no) = /i+(J, 0,no) — A, and by an identical argument for 
the hole excitations, that /i_(J, A,no) = /i_(J, 0,no) + A: the upper and lower 
pure system phase boundaries each simply translate inwards by A. This result is 
sketched in Fig. [5l The tip of the Mott lobe occurs at the point J — Jnj{^,no) 
where the two curves intersect, i.e., fi^{JM,0,no) — /i_ ( J^/, 0, no) = 2A, thereby 
replacing the critical singularity at Jc, A = by a simple slope discontinuity. 




Fig. 6. Schematic illustration of the droplet model of the Bose glass phase for the random site 
energy model (adapted from Ref. 27). Left: Exiting the Mott lobe at fixed J < Jm below its tip 
with increasing /i > fi+{J, A) [or decreasing /i < fJ.— {J, A)]. For arbitrarily small e = fi — fi^(J, A) 
there will be an (exponentially small) density of droplets with volume large enough that the energy 
gap to adding a particle is smaller than e. Right: Exiting the Mott lobe from its tip at fixed fi 
with increasing J > Jm- For arbitrarily small A J = J — Jjv/(A) there will similarly be arbitrarily 
large droplets whose local chemical potential magnitude is above the energy gap required to add a 
particle (previous droplets from the left illustration) , or below that required to remove a particle 
(additional droplets). 



4.2. Superfluid droplets and the Bose glass phase 

Now, as /i increases above /i+ (or decreases below /i-), these excitations begin to 
be populated with particles (or holes) — the system is compressible. As sketched in 
Fig. [6l the droplets now behave like uniform dilute superfluid regions just above or 
below the clean system phase boundary. As illustrated in somewhat more detail in 
Fig. [71 particles enter a given droplet at discrete values of increasing /i, with gaps 
controlled by the droplet volume and the bulk compressibility of the clean system. 
Although the droplets are superfluid, they are (exponentially in l/|/j, — /i±|) widely 

'More correctly, for larger J, one should replace Jd? in the second term of 5e by a quantity scaling 
with the background correlation length £,{,])■ 
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Droplet model schematic 




0.2 - |a^(J,A) = ~ 

o' ' ' ' ' 

5 10 15 20 25 

Droplet index i 

Fig. 7. Schematic illustration of spectrum of droplet excitations (adapted from Ref. I27j . Droplets 
are assumed to be well separated, and independent (see Fig.|6ll. The horizontal axis is an arbitrary 
droplet index. The vertical axis shows the critical chemical potential levels (in arbitrary units) at 
which a new particle is added for each droplet. The lower bound of the excitation spectrum lies at 
the Mott phase boundary /i^ (J, A) = fJ.+ {J, 0) — A (solid horizontal line). For any e = fi — /i+ > 
(dashed horizontal line), there will be a finite (but exponentially small) density of large droplets 
with local chemical potential lying sufficiently far above /i+ that one or more extra particles are 
added. The total number of particles added is given by counting the number of "occupied" levels 
below 11. For a given droplet, while the lowest excitation depends primarily on the local chemical 
potential, subsequent particles are added in sequence, with gaps scaling as 1/k+V where V is 
the droplet volume, and k+{J, fi+) is the bulk compressibility of the pure (superfluid) phase just 
above the Mott phase boundary. In an infinite system there are infinitely many droplets, implying 
a continuous distribution of excitation energies, and the bulk density will increase continuously 
with increasing e > 0. 

separated. Global superfluid long range order would require that the particles be 
able to tunnel coherently through the intervening Mott phase, producing a finite 
(but, at best, exponentially small) superfluid density. Such would indeed be the 
case, for example, for identical droplets arranged periodically in a homogeneous 
background. However, in the present case, the tunneling particle also experiences 
a residual random potential due to the Ui (suitably renormalized by many body 
effects in the Mott background), and the usual Anderson localization arguments 
imply that phase coherence cannot be maintained. 

A similar argument pertains if one exits the Mott lobe near its tip, by increasing 
J > Jm at fixed /i. In this case there are simultaneously present an exponentially 
dilute [in 1/(J— Jm)] set of both particle and hole superfiuid droplets (right panel 
of Fig. [5]). By the same Anderson localization argument, neither set of droplets can 
support global phase coherence. 

In summary, one concludes that a new compressible, insulting Bose glass phase 
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X 

Fig. 8. Schematic illustration of random rod-like structure generated by quenched disorder in 
quantum models (adapted from Ref. I27|l . The cylinders represent, for example, random regions of 
enhanced or suppressed hopping strength or site energies. 

completely surrounds the Mott lobe, and a direct MI-SF transition is impossible. 
However, as l/x — /i±| (or J — Jm) increases further, the superfluid droplets grow in 
size, and new ones are created (as /x overcomes larger and larger local Ui values). 
Eventually they percolate sufficiently that global phase coherence is stabilized, and 
a BG-SF transition takes place. This transition is presumably the ultimate zero 
temperature description of helium in Vycor. 



5. Path integral formulations and universality classes 

Given now a basic understanding of the underlying physics of the MI, BG and SF 
phases, I turn now to more formal treatments that allow one to clearly identify pos- 
sible universality classes, and their coarse-grained field theoretic descriptions. One 
begins with the path integral representation, based on the Trotter decomposition 
of the partition function of the Josephson function model Z = J D4)e^'~ ' with 
Lagrangian!23 

^ drl^Y. cos[^^i.r) - 0j(t)] + + i{u, - /i)]^ ) (5) 

in which /3 — \/kBT oo, K = 1/2V, and — oo < 0i(r) < oo are continuous 
classical phase fields. As usual, the quantum degrees of freedom give rise to an 
extra (imaginary time) dimension t in the effective classical model. In general this 
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extra dimension is highly anisotropic, and, except in special cases, does not simply 
lead to the same classical model in one higher dimension. It is immediately evident, 
for example, that the disordered quantities Jij, Ui are r-independent. The point like 
quenched disorder in the quantum model therefore leads to a picture of columnar 
or rod-like disorder in the classical model (Fig. [8]). 

In order to place ([5|) in the context of more familiar field theoretic models, we 
develop a coarse-grained, long wavelength continuum approximation to ([5]). Let 
gi0i(T) _^ 7/;(x, r), and relax the sharp condition jV'l = 1 on the field magnitude 
in the standard fashion by adding terms ^ro|-0p + iuo|'0|^ to C (which constrains 
IV'I ~ ro/uo). In addition the cosine term is approximated by a squared gradient 
term JjVV'p. The result is an effective ij}'^ Lagrangian, 

Cc[iP]= J^T I d'*x|i|VV'|2-iv*[a.-.9(x)]V + ^r(x)|V|' + JwolV^l'}. (6) 

Space and time have been rescaled so that the squared derivative terms have unit 
coefficient. The random is then subsumed into a random r(x) = rg + (5r(x), and 
—fj, + Ui has been subsumed into a random .9(x) = 170 + '55(x), with (dr) — {6g) — 0. 
Superfluidity corresponds to a nonzero anomalous average (ip) , and one may think 
of rg as the control parameter which takes one from the normal phase for large 
positive rg, to the superfluid phase for large negative rg. 

One may now consider the various possible classes of phase transition depending 
on the various terms, and their underlying symmetries, that one keeps in ([6]). Table 
[T] lists the relevant special cases. 

5.1. Particle-hole symmetric transitions 

The particle-hole symmetric model fi = 0, Ui = corresponds to ^(x) = 0. The 
Lagrangian Cq (clean system) or C2 (system with random hopping) is then purely 
real, and e~^"'^ may be interpreted in terms of a classical probability density for 
the two-component field (Re{^),lm{^)]. The squared-gradient interaction in both 
space and time implies that the model is the usual Landau-Ginzburg- Wilson rep- 
resentation of a ferromagnetic XY-model in {d + l)-dimensions. In the clean case, 
time and space are fully symmetric, and the transition (through the tips of the Mott 
lobes in Fig. [2]) is in the universality class of the usual classical {d+ l)-dimensional 
XY model. With random hopping, Sr{x.) 7^ 0, strong anisotropy is generated by the 
columnar disorder (see Fig. [8|) , and th e transitio n is in the universality class of the 
so-called classical random rod problem^^^^^ 

5.2. Particle-hole asymmetric transitions 

For nonzero fi or Ui, hence g{x), the Lagrangian is complex, and a classical proba- 
bilistic interpretation is no longer possible. The leading interaction is now a (purely 
imaginary) linear time derivative g'il)*drtp, which generates an an even stronger 
space-time anisotropy. For the clean system, go ^ but Sr{x) = Sg{x.) = (the 
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Pure PH-sym [{d + l)-diniensional XY model]: 

Co = J d'^X J dT {^\V^\^ + ^\drM' + 5^0|^P + jUom^} 

Pure PH-asym [d-dimensional dilute Bose gas] : 

PH-sym RR [{d + l)-dimensional classical random rod model]: 

PH-asym RR [{d + l)-dimensional incommensurate random rod model]: 
£3 = /d'^x/dr {i|VVP - '^ridr - go^i' + i[ro + 5ri^)M^ + luo\M^} 

Statistical PH-sym [commensurate dirty boson problem]: 
£4 = /d'^x/dr {i|VVP - '^r[dr - <5.9(x)]2V + i[ro + <5r(x)]|V'|2 + iuolVl^} 

Generic PH-asym [incommensurate dirty boson problem]: 
£5 = / d-^^; / dr {i|VV|2 - '^r[dr -go- 5gi^)]^i' + h[ro + Sr{^)M^ 

Table 1. Continuum 1/1* representation of models with various types of disorder and various degrees 
of particle-hole symmetry. The coefficients of IV^/ip and |9t'/'P have been normalized to i. The 
control parameters ro and go ^re analogous to Jq and //, respectively. Disorder in the hopping 
strengths is represented by 5r, while that in the site energies is represented by 5g. Both are 
independent of t, and in field theoretic treatments, are taken as quenched Gaussian random fields 
with zero mean and delta- function correlations characterized by variances and Ag, respectively. 
Disorder in the other parameters (including the unit gradient-squared coefficients) may also be 
introduced, but produces no new critical behavior. 



Lagrangian £1), this term generates critical behavior equivalent to the onset of 
superfluidity (described quantitatively by the Bogoliubov modeH at zero density in 
a dilute Bose gas0 

If 50 = and Sg{x.) ^ has an even distribution, then the corresponding 
Lagrangian £4 has a statistical particle-hole symmetry — the gip*dTip term is locally 
nonzero, but zero on average (the manifestation of this symmetry is more evident in 
the replicated, disorder averaged Lagrangian discussed in Sec. ^ . It will be argued 
below that the transition in this case is identical to that for the generic Lagrangian 
£5 in which all terms are nonzero, so that the statistical symmetry is restored at 
the critical point Q Intuitively, the expanding superfluid droplets in the Bose glass 

sSee, e.g., Ref. 1331 For a more modern view, see Ref. [14] 

''The usual coherent state functional integral formulation of the boson model JTJl corresponds to 
an absent term, and a unit coefficient ilj*dTip term, so that the model can never possess 

an exact particle-hole symmetry. Treating both particle-hole symmetric and asymmetric models 
requires the expanded space of Lagrangians corresponding to l6l l. 

'This idea has been used to explain the vanishing of the Hall conductivity at magnetic field-tuned 
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phase as the superfluid phase boundary is approached lead to a decouphng of the 
particle number from the lattice potential. The discrete droplet excitation spectrum 
pictured in Fig. [7l which replaces the individual lattice site occupation spectrum, 
becomes denser and denser, the statistics of the energy spectrum loses its up-down 
asymmetry, and there is less and less distinction between the statistics of particle 
and hole excitations as /i is varied. 

Less surprisingly, the particle- hole asymmetric random rod model, £3, in which 
go,Sr{x.) 7^ but 5g{:ic) = 0, is equivalent to the generic model. In the renormal- 
ization group sense, the combined effect of go and (5r(x) produces nonzero 5g{:ic). 
Intuitively, in the original boson model ([l]), random hopping produces regions of 
varying compressibility, and a uniform /i will produce a nonuniform density. As far 
as critical behavior goes, this is indistinguishable from density fluctuations produced 
a random site potential (Sec. [4] and Figs. [6l [7|). 

6. Quantum scaling theory 

In describing the Vycor data in Sec. [1] several zero temperature critical exponents 
were introduced. Here the general scaling theory of a QPT (which is not initially 
limited to the dirty boson problem) is summarized, these exponents are placed 
in a more general context, and scaling relations between them and more familiar 
exponents are derived.^^ 

Let S be the thermodynamic control parameter (e.g., J ~ Jc, fi — fic^ or — »'o,c), 
with critical point at 5 = 0. The divergence of the spatial and temporal correlation 
lengths, ^ w ^o|'5|^'^, ~ Co,t|^|~'^^ define critical exponents Vjiy-r- The dynamical 
exponent z = Vrjv quantifies the space-time anisotropy. In particular, the isotropic 
Lagrangian £0 must lead to z = 1, while for the remaining models one expects 
z 7^ 1. The singular part of the free energy scales in the form Fg w A|(5p^", which 
defines the analogue of a zero temperature "specific heat" exponent a. The quantum 
hyperscaling relation relation 2 — a = (d + z)!^ follows from the assumption that Fg, 
being an energy density, scale also inversely with the space-time correlation volume 

The critical behavior of the superfluid density, ps w ps^oli^r defines and exponent 
V. The superfluid density is a torsional modulus quantifying the stiffness of the 
order parameter, and is therefore derived from the free energy increment AFg = 
^PsVg (which must be part of Fg because ps vanishes in the insulating phase, and 
is therefore a singular quantity), where the superfluid velocity Vs = ftkg/m is 
generated by a long wavelength gradient — Vcf) in the order parameter phase 
(corresponding to a helical twist in the order parameter itself). Since kg is an 
inverse length, it must scale as and one obtains the quantum Josephson relation 
v = 2 — a — 2v= (d-|- 2: — 2)i/. If, as in the Vycor experiments, one uses the particle 
density difference as the control parameter, 5 = p — pc, then the exponent lu defined 



superconducting transitions: Refs. 1341 1351 
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in Fig. [T]is given simply hy uj ~ v. 

The temperature is similarly included via the temporal finite size scaling 
form Fs ~ A\6\'^~°'Y{(3/£^r), in which Y{y) is a universal scaling function with 
Y{oo) = 1. For finite P the argument is finite, and Y{y) must interpolate be- 
tween the zero and finite temperature critical behaviors. In particular, the finite 
temperature transition must occur at some value y — yc, which leads to the 
relation kBTc{S) « {yc£,T,o)^^\S\^'^ ■ Referring again to Fig. [U one therefore ob- 
tains ps{T = 0) ~ with 6 = {d + z — 2)/z. One obtains therefore the ratio 
Lu/9 — zv = Vr- 

Using the Vycor data experimental values uj = 1.7 ± 0.3 and 9 = 1.25 ± 0.2 
in d = 3, one obtains the reasonable result i^r = 1.4 ± 0.3. However, the implied 
relations z — 1/(0—1) and v = lo/Oz produce the remarkably unsatisfactory bounds 
2.2 < z < 20 and 0.07 < < 0.64. Unfortunately, more accurate bounds would 
require much higher quality data, but even with the intervention of 25 years, the 
original experiments - would be extremely difficult to improve upon. 

6.1. Compressibility and z vs. d 

If one uses 5 ^ /i — /ic, the compressibility k — dn/dfi ~ —d^F/dS^ has a singular 
part Ks ^ controlled by the zero temperature specific heat exponent (this 

is true more generally so long as variation in S leads to variation in the density). 
However, the total compressibility also has an interpretation as a temporal superfluid 
density: the free energy increment due to temporal phase twists LUg — drip is given by 
AFr = ■^Kuj'^. This identification follows from (O since the imposition of the phase 
twist ujs is equivalent to the shift fj, ^ ^ + icug, and it follows that d'^ AFr / duj'^ = 
—d'^F/dfi^ = K. By analogy to ps, if one assumes that AFr is part of Fs, and 
proposes that Ug scales as then one obtains k ~ Kolf^T^ with = (d — z)i/. 
The physical fact that k is finite and nonzero through the transition (with Kg 
yielding a small correction since one expects a < 0) then leads to the proposed 
scaling relation z = d^ 

However, this argument lies on very shaky ground because it is difficult to 
reconcile Kg and k obeying separate critical scaling relations. In fact, the difference 
K — Kg, being finite through the BG-SF transition, is dominated by the analytic 
background, and AFr should be likewise.'SS] Moreover, since the shift ^ + iujg 
simply leads to a small adjustment of a parameter that is already nonzero, the 
scaling assumption ujg ~ is inappropriate. Only if p, = 0, Ui = 0, so that 
ajg breaks particle-hole symmetry, should one expect such a scaling^ Indeed, for 
the random rod problem, t\i — t\i g IS entirely singular since, like pg, it vanishes 
identically in the RRG phase. Consistently, one finds Vr > 0, hence z < d, at that 

jPor ps, the spatial phase twist corresponds to the shift — » (V+iks)^, introducing a momentum 
term ik^ • ?/)* Vi/' which breaks a spatial inversion symmetry. It is for this reason that one expects 
the critical scaling of kg with 
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transitionPSEHSSl Consistently as well, the control parameter 5 = J — Jc in that 
case does not couple to the density, and k is no longer related to the exponent a. 

In conclusion, at both the RR-SF and BG-SF transitions z is expected to remain 
an independent exponent, undetermined by any simple scaling relation. Recent high 
resolution quantum Monte Carlo simula tion s in d = 2 support this view, finding 
z = 1.40 ± 0.02 at the BG-SF transition.ISS 

6.2. Universal critical sheet conductance in d — 2 

The conductivity is related to the superfluid density via the zero frequency limit 
of cr = Ps/iu). If one assumes that w, being an inverse time, scales as (this is 
reasonable here, because w is a real dynamic frequency, not an imposed phase twist, 
and a is indeed singular), one obtains a w ctqIiJI^ with E = 2—a~{z+2)v = {d—2)v. 
This has the remarkable implication that, very generally, in c? = 2 the conductivity, 
while vanishing in the insulating phase, and diverging in the superfluid phase, has 
a finite constant value CTc right at the critical point. Experiments on disordered 
superconducting thin films, wher e i5 is the film thickness or an applied magnetic 
field, provide some support for thiJ^HESHO] ^ggg g^jg^ Ref. 1411 and references therein). 
Moreover, quantum hyperuniversality^^ which predicts not only that the critical 
exponents for Fg and {£,t£,'^')^^ coincide, but also that the amplitude combination 
Fg^r^'^ ~ ^?T,oCo s-lso be a universal number, predicts that ctc — A^r,o£,o coincides 
with this amplitude, and therefore should be MnwersaZ0 This has much les s exper- 
imental support — a large range of estimated ctc values are found .I^^^^IHTI ^j^g 
other hand, critical amplitude ratios are far more difficult to measure accurately 
than critical exponents, and higher quality data would again be desirable. 

7. Renormalization group approaches 

The dirty boson problem has so far discussed from a phenomenological point of 
view, emphasizing elementary excitations, symmetry principles, and results that 
follow from general scaling relations. This review is concluded with a brief survey 
of renormalization group approaches to quantitative evaluation of exponents and 
other quantities. 

7.1. One- dimensional models 

There is one limit where an exact solutions exist, namely d = i l 44 | 45 | 7 | 28 | 29 | rpj^^ 
model maps to a classical 2D fluctuating sine-Gordon-type interface roughening 

''The scaling theory is general, and does not require any identification between the dirty super- 
conducting and superfluid systems. However, in modeling thin films close to the transition, it is 
typically argued that Cooper pairs may indeed by treated as bosons, hence that the transition lies 
in the same universality class as the dirty boson problem, and Cc may be computed within this 
model. See, e.g., Ref. 1431 
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niodel!23 

/:sG[h]- JdxJ dTl^^K{x){drhf + ^V{d^h)^-m,{x)d^h-2yoCOs{2Trh)Y{7) 

where h{x,T) is the interface height at a given space-time point The hopping dis- 
order leads to disorder in the interface tension coefficient K{x), while the chemical 
potential and site disorder leads to the disordered tilt potential m{x) = mQ + Sm(x). 
The cosine term prefers integer values of h, a reflection of the original integer boson 
site occupancy. At sufficiently large yo this term pins the interface at a fixed integer 
height, with small fluctuations about it. This is the Mott insulating phase. For the 
ID random rod problem, m{x) = 0, at small yo the fluctuations are able to unpin 
the interface, leading to logarithmically divergent large scale fluctuations: this cor- 
responds to the superfluid phase. The disorder in K turns out to be irrelevant at 
this transition, so both the clean and random rod problems are described by the 
usual classical 2D Kosterlitz-Thouless theory, with, in particular z = 1 and critical 
correlation decay exponent t] = j (see, e.g., Ref . l47| . 

On the other hand, nonzero ni{x) couples directly to the interface slope, and 
can have a much stronger effect. At a sufficiently large average value, mo, there 
is a transition to an interface with an average constant spatial tilt. The value of 
mo at the onset of this tilt corresponds to the Mott excitation gap. For sufficiently 
small K, V the fluctuations about the mean tilt completely wash out the cosine 
term and are again logarithmically divergent on large scales: this is again the su- 
perfluid phase. However, at intermediate values there is enough residual effect of 
the cosine term to constrain the interface to finite fluctuations about this tilted 
state — this corresponds to the Bose glass phase. The BG-SF transition is described 
by a modified Kosterli tz-Thouless theorj^j which produces different universal crit- 
ical properties Yot example, although z = 1 still, one now finds 77 = i. 
Moreover, near the transition one may absorb mo into a shift h h — (rno/V)x, 
so that only Sm enters. This explicitly demonstrates the asymptotic restoration of 
statistical particle-hole symmetry in the ID model. 

7.2. Epsilon- expansions in higher dimensions 

We now turn to epsilon-expansion type approaches in higher dimensions that have 
proven so useful in quantifying classical critical phenomena. Renormalization group 

'An extension to d = I + e dimensions has also been proposed, based on inserting the engineering 
dimensions appropriate to higher d into the one-dimensional RG flows: see Ref. 46 However, in 
the absence of a form for the surface roughening type Hamiltonian generaUzing ^ for noninteger 
d > 1 (the particle- vortex duality transformations used to obtain such models requires integer d), 
there is presently no rigorous support for this approach. 

™In the replicated model, the tilt disorder replaces the usual cos{2iTh) term by a replica interaction 
term of the form cos{[27r[/ia(a;, t) — /i^(a;, t')]}. This term generates renormalization group flows 
similar to those of classical Kosterlitz-Thouless (Ref. 147)1 but with sufficiently different geometry 
that different universal critical properties are produced. 
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approaches are generally based on "replicated" Lagrangians, in which the disorder 
has been integrated out at the expense of introducing p copies, ipa, a = 1, 2, . . . ,p 
of the original field that all interact with the same quenched disorder Sr,Sg, but 
are otherwise independent, with the formal limit p —> taken at the end. Using the 
standard delta-correlated Gaussian disorder model, with variances Ar and Ag, the 
replicated -0^ Lagrangian takes the form C^'' — Cq\ + Cq\ withj^^ESl 



E f dr f dr' [ d''x{Ar\M^,T)\'\M^,T')\^ 

+ A3[Ca.Va-5o|^a|'](x,T)[V^^5,V,3-5o|V',3|'](x,r')}. (8) 

The disorder average leads to quadratic interactions between the replicas that are of 
infinite range in t, a reflection of its columnar character, and this is what generates 
new critical behavior. The random rod model corresponds to 50 = 0, Ag = 0, so 
that the Lagrangian is separately invariant under reversal of t and r'. The case of 
statistical particle-hole symmetry corresponds to 50 = 0, Ag > 0, which leads to a 
symmetry only under simultaneous reversal of t, r'. 

Now, imagine beginning with the random rod model, and perturbing it in var- 
ious ways. Nonzero Ag leads, in a heuristic Fourier space notation, to a Aga;^'0^ 
correction to the Ari/j^ term already present. Under naive power counting, since 
critical fluctuations are dominated by large scales and long times, increasing powers 
of frequency lead to less important terms, and one would expect the former term 
to be strongly irrelevant relative to the latter, i.e., that symmetric site disorder 
should not destabilize the RR critical behavior. On the other hand, go introduces 
a igoujip'^ perturbation, which by power counting dominates the tu^il''^ term, and is 
expected to lead to new critical behavior. This is certainly conflrmed in the clean 
model, where it generates the crossover from {d + l)-dimensional XY behavior to 
weakly interacting Bose gas behavior (Sec. I4.1| l. However, the notion of asymptotic 
restoration of particle-hole symmetry requires that, in the disordered case, this term 
ultimately be subdominant to the Ag term. 

These paradoxes find their resolution in the fact that naive power counting is 
exact only at Gaussian critical points, and approximately valid only when interac- 
tions (terms of order higher than ■0^) play a small role. If their role is large, then 
different terms can indeed exchange dominance. In order to explore such notions 
quantitatively, it is convenient to have a model parameter that, when varied, inter- 
polates between the near-Gaussian and strongly non-Gaussian limits. In classical 
critical phenomena, this parameter is e = 4 — d: the standard epsilon expansion is 
based on the fact that the critical behavior is Gaussian for d > 4, and nearly so 
for small e > 0. The corresponding approach for problems with columnar disorder 
requires not only that d be close to four, but also that the dimension Cr of the 
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go 




Fig. 9. Proposed behavior of the random rod (RR), commensurate (statistically particle- hole 
symmetric) dirty boson (C), and incommensurate (fully particle- hole asymmetric) dirty boson 
(IC) fixed points as functions of e-r- (adapted from Refs.l27l l48b . The axis has been suppressed. 
Here G and Gg are commensurate and incommensurate Gaussian fixed points. The commensurate 
fixed point bifurcates away from the random rod fixed point at ~ At er = eri ~ ?, C 



and IC merge and at e-r 
problem. 



1, C is the stable fixed point that describes the physical dirty boson 



columns be small — essentially this is an expansion about the point-disorder limit. 
With D = d + er being the total dimensionality, and e = 4 — D, one may develop 
renormalization group recursion relations for the random rod problem which ha ve a 
perturbatively accessible critical fixed point when e, e-r are both small.lSSEIlSH xhe 
physical problem in d = 3 is recovered for Et = 1 and e = 0. 

This double epsilon-expansion is very poorly behaved — with coefficients having 
very strong dependence on these parameters, and higher order terms dominating 
lower order ones unless e, e^- are extremely small. Nevertheless, as shown in Fig. 
ini a very compelling picture emerges. For sufficiently small e^- the naive power- 
counting expectations are borne out: the random rod fixed point (RR) (which be- 
comes the usual Gaussian fixed point G at e = e,- = 0) is stable against statistically 
particle- hole symmetric perturbations {go = 0, but Ag > 0). However, it is strongly 
unstable to nonzero go^ which generates a crossover to a new, fully stable, particle- 
hole asymmetric, or incommensurate (IC), critical fixed point (which becomes a 
Bose Gaussian fixed point Gb at e = e,- = 0). 

However, for Cr > (with = ^ to leading order in e^), nonlinear terms in 
the recursion relations dominate the linear ones (the latter reflecting naive power 
counting), and the RR flxed point becomes unstable to A^, and a new statistically 
particle- hole symmetric, or commensurate (C), fixed point bifurcates away from 
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RR. For not too large Cr both C and RR remain unstable to go 7 and IC hence 
remains the globally stable critical fixed point. However, for e^- > e^i (with eri — | 
to leading order) the incommensurate fixed point intersects the go = plane, and 
merges with C. The latter is now completely stable, and is proposed to correspond 
to the dirty boson fixed point at Cr = 1 • This provides a detailed scenario by which 
statistical particle-hole symmetry is restored. However, the fact that e^, though 
small, is a finite number means that the fixed point C is not perturbatively accessible 
in the usual sense, and one should treat the detailed estimates in the vicinity of 
this fixed point with caution. The value e-ri is even larger, and therefore even more 
uncertain, and extrapolation of these results to Cr = 1 should be treated as, at best, 
qualitative estimates. The general scenario proposed, however, seems very natural 
and illuminating. 



8. Summary and conclusions 

Over the past twenty years, advances in understanding the rich physics of the 
phases and phase transitions in disordered boson systems has proceeded fruitfully 
along a number of d ifferent fronts, including descriptions of excitati ons w ithin the 
different phases scaling p henomenology for the critical expo nent j^^*^^! and uni- 
versal amplitudes,!^^^^^^^!!!! exact solutions in one dimension,''^^''^^- ^^ ! ^^ ! ^^ ! nu- 
merical simulations i n one and two dimensions ,-1^ '^^ ^S! weak hop- 
ping expansion work,f23l24j ^^^^ approxir nate r enormalization group epsilon ex- 
pansion approaches in higher dimensions .^^ZEHl j have tried here to give a flavor 
for all of these, emphasizing qualitative features, such as phase diagrams and 
the impor tance of particle -hole symmetry, ' ° and making connections to 
experimentl ^ l ^^ l ^^ l '^^ l '^^ l ^ ^ where possible. Perhaps the most interesting future de- 
velopment would be quantitative estimates of the BG-SF critical behavior in three 
dimensions, which appears to be quantitatively accessible only through numerical 
simulation. This would allow direct comparisons with the Vycor data (at least for 
those exponents that appear to be reasonably constrained by the data; of course, 
better experimental helium data, perhaps using other disordered substrates, would 
be hugely beneficial as well). However, this may have to await future generations 
of computational capability, as the required system sizes are currently well out of 
reach. 
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